Isometric and metamorphic operations on the space of local 

fundamental measures 



Matthias Schmidt^' ^ 

^ Theoretische Physik II, Physikalisches Institut, 
Universitdt Bayreuth, D- 9 5440 Bayreuth, Germany 
^H. H. Wills Physics Laboratory, University of Bristol, 
Royal Fort, Tyndall Avenue, Bristol BS8 ITL, United Kingdom 
(Dated: 15 November 2010, revised: 23 December 2010, to appear in Molecular Physics) 

Abstract 

We consider symmetry operations on the four-dimensional vector space that is spanned by the 
local versions of the Minkowski functionals (or fundamental measures): volume, surface, integral 
mean curvature, and Euler characteristic, of an underlying three-dimensional geometry. A bilinear 

combination of the measures is used as a (pseudo) metric with + H signature, represented by a 

4x4 matrix with unit entries on the counter diagonal. Six different types of linear automorphisms 
are shown to leave the metric invariant. Their generators form a Lie algebra that can be grouped 
into two mutually commuting triples with non-trivial structure constants. We supplement these 
six isometric operations by further ten transformations that have a metamorphic (altering) effect 
on the underlying geometry. When grouped together, four different linear combinations of the 
metamorphic generators form a previously obtained third-rank tensor. This is shown to describe 
four different types of mutually commuting "shifting" operations in fundamental measure space. 
The relevance for fundamental measures density functional theory is discussed briefly. 
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I. INTRODUCTION 



Applying classical density functional theory (DFT) requires to have an approximation for 
the Helmholtz free energy as a functional of the one-body density distribution(s) jl, 2|. For 
;he case of additive hard sphere mixtures, Rosenfeld's fundamental measures theory (FMT) 
sl is an approximate DFT that unified several earlier liquid state theories, including the 
Percus-Yevick integral equation theory and scaled-particle theory, and encapsulates their 
results in a free energy functional. Several recent reviews give a detailed account of FMT 
and some of its extensions and modifications {4-^. The theory was used to address a broad 
variety of interesting equilibrum phenomena, ranging from freezing to capillary behaviour 
of liquids. When compared to computer simulation data, theoretical results for e.g. density 
profiles and interfacial tension were typically found to be very reliable. FMT rests on building 
weighted densities via convolution with the bare density profile (s). The microscopic density 
profile Pi{r) of species i gives the mean number of particles of species i in an infinitesimal 
volume element at given position r and hence carries dimensions of (length)"^. The weighted 
densities in FMT are smoothed versions of these "real" density distributions. In Kierlik and 
Rosinberg's (KR) version of FMT j^, there are four scalar weight functions for each 
hard sphere species. Rosenfeld's original approach that involves additional weight functions 
was shown to be equivalent to the KR version js], and was later carried much further by 
Tarazona ^ and Cuesta et al [lO|. FMT has intimate connections to methods from integral 



geometry [Uj via the Gauss-Bonnet theorem [l2l |. 

The weight functions in FMT are quantities with dimension of negative integer powers 
of length, ranging from (length)"^ to (length)". A linear combination of pairs of weight 
functions that are convolved with each other is used to express the Mayer bond fij{r), as 
a function of distance r. Recall that for hard sphere mixtures, the Mayer bond equals 
fij{r) = — 1 for distances r < Ri + Rj, i.e. when the two spheres with radii Ri and Ri over- 
lap, and it vanishes otherwise. Here the subscripts i,j label the different species. Originally 
proposed for hard spheres, this framework was sufficient to to derive FMTs for models such 

n n 

as the Asakura-Oosawa colloid-polymer mixture [13| and the Widom-Rowlinson model |14| . 
However, the treatment of binary non-additive hard sphere mixtures required significant 
modification of the mathematical structure of FMT In particular, further weight func- 
tions were introduced in order to correctly model the deviation of the hard core interaction 
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range between species i and j from the sum of their radii, Ri + Rj. The fact that this 
deviation is non- vanishing is the defining feature of non-additive hard sphere mixtures. The 
additional weight (or kernel) functions possess dimensionalities up to (length)"^. They can 



be grouped in a double-indexed tensorial form [15| and were shown to possess a remark- 
able group structure [l6|. In very recent work, the FMT for non- additive hard spheres was 



applied successfully to bulk structure [l7|, ll8| and to interfacial phenomena [18 1. 



Several features of the mathematics that underlies the FMT weight functions have 
emerged Q Q: i) The four different position-dependent fundamental measures (in the 
KR formulation) can be viewed as elements of an abstract four- dimensional vector space, ii) 
Based on dimensional analysis, a (pseudo) metric can be defined, which can be represented 
by a 4x4 matrix with unit entries on the counter-diagonal. All other entries in this matrix 

vanish. The metric has a + H signature, hence it differs both from that of Minkowski 

spacetime in special relativity (H ) and from that of four-dimensional Euclidian space 

(+ + ++). iii) Operations that are common in linear algebra, i.e. matrix multiplication 
and more general contraction of tensor indices possess meaningful interpretation, see e.g. 
the shifting transform described in Ref. [3]. Here all product operations are carried out in 
Fourier space and hence correspond to convolutions in real space. 

In the present paper we explore the mathematical structure further by focusing on sym- 
metry operations that leave the fundamental measure metric invariant. Our motivation 
comes from the fact that careful analysis of the symmetries is central to exploiting the prop- 
erties of any (abstract) space. Typically, this tasks requires the identification of the linear 
automorphisms that leave the metric invariant. Recall that an automorphism is a bijective 
function that maps a space onto itself (i.e. both function value and argument are element 
of the same space). Much structure can be revealed by considering infinitesimal versions, or 
generators, of the transformations. In a Lie algebra the commutator of any pair of generators 
can be represented as a linear combination of again the same generators. The coefficients of 
the linear combinations form the structure constants of the algebra. 

In Euclidian space, the symmetry operations that leave the metric invariant are orthogo- 
nal transformations, or rotations. These possess three (six) independent generators in three 
(four) spatial dimensions. For the case of Minkowski spacetime with three spatial and one 
time-like dimension, there are three spatial rotations and three Lorentz transformations, or 
boosts, the latter coupling time and one of the spatial dimensions. The number of degrees 
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of freedom, and hence the dimensionahty of the group of isometries, is independent of the 
signature of the metric. However, the algebraic structure, as expressed by commutator re- 
lations between the respective (infinitesimal) generators of the transforms, differs for both 
cases. For the case of spacetime, the resulting mathematical structure is the Lorentz group. 
(One refers to the Poincare group when four translations in the different spacetime direc- 
tions are added.) Here we present in detail a similar analysis for the space of Minkowski 
functionals jll|. We describe four boosts and two rotations that leave the metric invariant. 
These are complemented by further ten operations that change the metric and that we refer 



to as metamorphic operations. We show that the spherical shifting operation of Ref. [16| is 
readily generalized to four different types of shifting, and that the corresponding generators 
can be expressed as linear combinations of the metamorphic generators. 

The paper is organized as follows. In Sec. Ull the theory is laid out, including the de- 
scription of inner boosts and inner rotations as isometric transformations (Sec. Ill Bp . of 
metamorphic operations (Sec. Ill C[) . and the relationship of Jeffrey's third-rank tensor jicl 
to the latter (Sec. IIIDI) . Concluding remarks are given in Sec. IIIII 



II. TRANSFORMING THE FUNDAMENTAL MEASURES 



A. Metric and inner scalar product 

We consider a four- dimensional real vector space with elements u = {uq, ui, U2, u^), where 
u depends on the three-dimensional argument q in Fourier space. The dependence on three- 
dimensional position r is then obtained by inverse Fourier transform, (27r)~^ J dqe^^ '^u. The 
vector components u^, with index u = 0,1, 2, 3, are dimensional objects: u^, possesses the 
dimension (length)''. Hence U3 is a measure of volume, U2 of surface, ui of mean curvature, 
and no of Gaussian curvature. We let the take on arbitrary (real) values, and hence 
restrict ourselves not to cases where the measures describe an underlying geometrical body. 
The interpretation of the u^, in terms of geometric measures is only intended to guide the 
intuition, the mathematics that we present in the following is based on formal arguments. 
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We use the metric represented by the matrix 

/O l\ 

10 

M = , (1) 

10 

\l 0/ 

hence a measure of squared "length" of a vector u is given by • M • u = 2(uoU3 + U1U2), where 
the superscript t indicates matrix transposition, and the dot indicates matrix multiphcation. 
The scalar product between two vectors u and v is • M • v = U0V3 + U1V2 + U2V1 + U3V0. 
Clearly, this is symmetric upon interchange of the vectors, i.e. ■ M ■ v = ■ M • u. The 

eigenvalues of M are —1 and 1, both are doubly degenerate; hence M possesses (+ H ) 

signature. As M is not positive definite (i.e. not all of its eigenvalues are positive), it can 
yield negative squared distances and hence constitutes not a metric in the strict sense, but 
one refers to a pseudo metric. While it is enirely possible to discriminate between covariant 
and contravariant vectors and correspondingly introduce lower and upper indices, which can 
be interchanged by application of the metric, we will not do so in the following. The present 
paper is primarily concerned with second-rank tensors, and we find the (index-free) matrix 
notation to be simpler, and will primarily rely on this in what follows. 

The hard sphere weight functions of FMT can serve as an example. These are functions 
of the squared wave number and the radius it! of the hard spheres. The Fourier space 
expressions of the KR version of the weight functions are wq — c-\- {qRs/2), wi = {qRc + 
s)/{2q), W2 = AirRs/q, = An{s — qRc)/q^, where s = sin(gi?) and c = cos{qR). The 
real space expression that corresponds, via inverse Fourier transform, to w-i, is a unit step 
function with range of R, i.e. Q{R— |r|), where 6(-) is the Heaviside (step) function. Within 
our framework we view the Wi, as the four components of a vector w. By straightforward 
explicit algebra one can show that • M • w yields the Fourier transform of a unit step 
function with range 2R, i.e. the expression for given above, but with R replaced by 
2R. Explicitly this is 2{woW3 + W1W2) = A7i[sm{2qR) — 2qRcos{2qR)]/q^. The significance 
in statistical physics stems from the fact this is (up to a trivial minus sign) the Fourier 
transform of the negative Mayer function of the pair potential of hard spheres of radius R. 
For a mixture, the additional species possesses weight functions v^, of range R', given by the 
above expressions for Wu, but with it! being replaced by R'. It is straightforward to verify 
that V M • w = 47r[sin(g(i? + R')) - q{R + R') cos{q{R + R'))]/q^, which again is the above 
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expression for the unit step, W3, but with R replaced by the sum of the radii, R + R'. These 
identities constitute one of the central building blocks of KR's formulation of FMT. The 



generalization to non-additive mixtures [15| amounts to introducting 4x4 matrices that 
change the range of the weight functions w,^. This shifting operation is discussed in detail 
in Ref. 16|]. Below in Sec. IIIDI we give further three such "internal" shifting operations. 
We emphasize that all transformations that are considered here are of internal nature, i.e. 
act on the four-dimensional space of fundamental measures, as opposed to e.g. translations 
and rotations of the underlying three-dimensional Euclidian space, which we do not consider 
here. 

The central aim of this paper is to formulate linear automorphisms that leave the met- 
ric ([T]) invariant. We refere to such operations as isometrics on the space of fundamental 
measures. Hence one has to identify 4x4 transformation matrices A that obey 

A* . M ■ A = A ■ M ■ A* = M, (2) 

which implies that a vector u and its transform A • u possess the same squared modulus. 
This can be seen from (A ■ u)* • M • (A ■ u) = • A* ■ M ■ A ■ u = ■ M ■ u, where the last equality 
follows from (|2]). An alternative is obtained by multiplying (|2]) from the right by the inverse 
A^^, and from the left by M, and observing that = 1, where 1 is the 4x4 unit matrix. 
Hence 

M-A*-M = A-^ (3) 

Note that this differs from the condition for orthogonal matrices. A* = A^^. While transpo- 
sition can be viewed as mirroring the matrix elements on the diagonal, the operation on the 
left hand side of corresponds to mirroring the matrix elements on the counter diagonal. 



B. Inner rotations and boosts as isometries 



Let us formulate the linear isometries, i.e. the automorphisms A that obey ([2]), by choosing 
appropriate generators X for each different type of transform, where X is a 4 x 4 matrix. 
The transformation matrices A are then obtained by (matrix) exponentiation. In order to 
see this, consider that the expression 1 + Xdr can be viewed as an infinitesimal transform 
of differential magnitude dr. A transform by a finite amount r can then obtained in the 
continuum limit of A^-fold application of the infinitesimal transform, where each step is taken 
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to be of magnitude t/N. This amounts to limjv->.oo(l + rX/N)^ — exp(TX) = A, where 
the result depends parametrically on r and the form of A is specific to that of X. Here 
the exponential of a matrix is defined by its power series exp(rX) = X]m=o(^'"/"^')^'"- 
the following, we allow the transformation parameter r to be dimensional, i.e. to carry a 
non- vanishing power of length scale. 

In order to allow for meaningful matrix multiplication (as is necessary for matrix expo- 
nentiation) the generators need to possess matrix components with suitable dimensionalities. 
This implies that the product of the transformation parameter and a matrix entry, rX^j,, 
where n enumerates the rows and u enumerates the columns, with both indices running 
from to 3, must be of unit (Icngth)^^*^. Taking matrix powers then preserves the ordering 
of dimensions, i.e. the yUZ/-component of the m-th matrix power, (t^X"*)^!,, has the same 
dimensionality as rX^,^ itself. Hence we can exponentiate the generators and obtain finite 
transforms. Besides letting r be a dimensional object, in the following the only further 
dependence on length scale shall be via q^, the squared argument in Fourier space. This 
corresponds to the (negative) Laplacian in the corresponding real three-dimensional space. 

From general arguments for four- dimensional spaces, we expect the isometry group to be 
six-dimensional, i.e. to possess six linearly independent generators, cf. the cases of Euclidian 
space and Minkowski spacetime of special relativity mentioned above. Given a set of such 
generators, {X^}, ennumerated by index a, a general transform is obtained as exp(^^ t^Xq), 
where is the magnitude of the ct-th transform. In principle the different contributions 
to the the total transform can be disentangled via the Baker-Campbell-Hausdorff formula. 
This requires knowledge of the algebraic group structure, which is encoded in commutator 
relations between the different generators, as laid out below. 

Here we discriminate between four generators for boosts, Bq,, and two generators for ro- 
tations. Da- The subscript indicates the dimensionality; the /lu-element of a given generator 
matrix possesses units of (length) As laid out above, all elements along a given di- 
agonal possess the same dimensionality; the dimensionality then decreases (increases) by 
one power of length scale when moving up (down) to the next diagonal. We call boosts 
those generators that satisfy ■ Bq, = Generator of rotations are those that satisfy 

Da ■ = — The generators are not unique; one can always build linear combinations 
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to obtain a different formulation. Here we choose the following set of generators. 

/O -q^ 0^ 

10 

\o -1 y 

/O -g2 0\ 

10 
g2 

\0 -10/ 

We have grouped the generators into two families, each consisting of two boosts and one 
rotation. The first one consists of Bq, B2, D2 and is given in (jl]), the second one consists of 
Bo', Bi, Di and is given in ([5]). Both families form closed Lie algebras, constituted by the 
commutator relations 



(5) 



[Bo, B2] = 2D2, [Bo, D2] = 2B2, [B2, D2] = 2g^Bo, 
[Bo',Bi] = 2Di, [Bo',Di] =2Bi, [Bi, Di] = 2g2Bo., 



(6) 
(7) 



where the commutator between two matrices X and Y is defined as [X, Y] = X ■ Y — Y ■ X. 
Members of different families commute; these are pairs of boosts: [Bq, Bq'] = [Bq, Bi] = 
[Bo',B2] = [Bi, B2] = 0, the (only) pair of rotations: [Di, D2] = 0, and the four mixed 
pairs of a rotation and a boost: [Bo, Di] = [Bi,D2] = [B2, Di] = [Bq/, D2] = 0. Tab. [J 
gives an overview of the group structure in table format. All relationships can be obtained 
by straightforward matrix algebra. We give a full multiplication table in Tab. [TTl anti- 
commutator relations are included for completeness. Note that in each family already the 
bare products (not commutators) give the result of the commutators up to a factor of 2. As 
a consequence, the anti-commutators within each sub-algebra vanish, see Tab. [H Altough 
commutators between members of different sub-algebras vanish, their plain products do not, 
cf. Tab.Hn The nine matrices that result from the products (as referred to in the off-diagonal 
blocks in Tab. [Tll) will be used below in order to define further, metamorphic, operations on 
the fundamental measures. 

It is now straightforward to calculate finite transforms via exponentiation of the respective 
generators multiplied by its transformation parameter. Recall that the latter is a dimensional 
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TABLE I: Table for commutator relationships [X, Y]/2 for the generators of boosts, 6,^, and rota- 
tions, Djy. X denotes a matrix of the leftmost column, Y one of the top row. 
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TABLE IL Left: Multiplication table X ■ Y for products of the generators of boosts and rotations. 
Right: Table of anti-commutator relationships {X, Y}/2 for the generators of boosts and rotations. 
In both tables X denotes a matrix of the leftmost column, Y one of the top row. 

object, and that the most general finite transform is given by exp(^j^ Ti^Xj,), where Ti, possess 
dimensions of (length) Here we give only the results for the case where all parameters bar 
one vanish. These are the following expressions for finite transformations corresponding 
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(9) 



(10) 
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g^ sinh(r2g^ 


g" sm(^r2g' 


cos(r2g^) / 

When applied to a vector u, (|8]) describes a multiplication of the components mq and u\ by 
e^°, and division of M2 and U3 by the same constant. Trivially, the (pseudo) squared modulus 
2(mi'U2 + ■Uo'^s) is left unchanged. Eq. OH]) is reminiscent of a hyperbolic rotation, and ( ITU]) of 
an ordinary rotation. Note the difference in occurrence of the minus signs in ^ and f[TII|) . 
For the generators ([5]) we obtain the following finite transforms: 
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cosh(rig) gsinh(rig) 
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(12) 



(13) 
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Again (fTTj) induces a straightforward scaling of vector components, ( !T2|) is a hyperbolic 
rotation and is f|T3l) an ordinary rotation. Recall that hyperbolic rotation can be viewed as 
Lorentz transforms (and vice versa). 

As a summary, we have identified six real matrices Bo, Bq', Bi, B2, Di, and D2, that posses 
the algebraic structure shown in Tab. [ITl The general (real) linear group, i.e. that of all real 
4x4 matrices, is 4^ = 16 dimensional. Besides the unit matrix, this leaves nine matrices to 
be considered. In the following we will use the matrices obtained as products of two isometric 
generators, cf. Tab. [Tll We find it interesting to investigate their action, when viewed as 
infinitesimal transformations, on the space of fundamental measures. Clearly, they cannot 
generate isometrics - we have exhausted these already. Hence we expect that the metric 
will not be conserved under the application of these further transformations, and we will 
henceforth refer to these transformations as metamorphic, as they change the underlying 
geometry in a fundamental way. 

The difference between automorphism and metamorphisms is reflected in the symmetry 
properties of their generators. The isometric generators (HI) and (jS]) are anti-symmetric with 
respect to mirroring on the counterdiagonal, i.e. each generator X satisfies 



This can be seen by inserting the infinitesimal versions A = 1 + Xdr and A^^ = 1 — Xdr 
into ([3]). Note that the symmetry f|T^ leaves 6 parameters free, which is consistent with 
the dimensionality of the corresponding group of transformations (and hence the number of 
generators). Correspondingly, metamorphic generators are symmetric under mirroring on 
the counter-diagonal, i.e. they satisfy 



as we will see in the following. Note that the symmetry ( IT5|) leaves 10 parameters undeter- 
mined. 

C. Metamorphic transformations 

We start by giving the explicit expressions for the matrices that we choose as generators 
of the metamorphic operations. As above, the index u of a given generator indicates its 



M ■ X* • M = -X. 



(14) 



M ■ X* ■ M = X, 



(15) 
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dimensionality. Explicit expressions for the nine different generators are as follows. 
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Here we have grouped the nine generatores into three Abelian subgroups, given in (fT6|) . 
f lT7|) . and f|T8|) . respectively. Any pair of matrices from of one of these subgroups satisfies 
[X^,Yj,] = 0. In general, the commutator between matrices from different subgroups is 
(up to a minus sign) a multiple of g^ times an isometric generator. Some of these pairs 
commute. The complete algebra of commutator relationships between the metamorphic 
generators is summarized in Tab. IIIII Remarkably, the commutator between any two pairs 
of these matrices either vanishes or it is a multiple of one of the isometric generators of 
Sec. IIIBI Inevitably, some of the quite compact structure of the previous subsection is lost, 
due to the sheer number of possible pairs. Nevertheless, note that indeed members of the 
same triplet {Fi, F2, F3}, {Hi, H2, F3/} and {Pq, P3, P3'} commute with each other. We defer 
multiplication and anti-commutator tables to the appendix. Clearly the nine generators are 
not unique. In the following section we will relate a previously obtained third-rank tensor 
to a linear combination of these generators. Before doing so we give explicit expressions for 
the finite metamorphic transformations. 
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TABLE III: Commutator relations [X, Y]/2 for the generators of metamorphic and isometric oper- 
ations. X denotes a matrix of the leftmost column, Y one of the top row. The upper block of nine 

rows give the commutator between pairs of metamorphic generators. The lower block with 6 rows 
give the commutators betweenone isometric and one metamorphic generator. 
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For the first triplet of generators, these are 
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\—q ^sinh(T3g^) 











The second group of finite metamorphic operations is 

^ cosh(TiQ') — gsinh(Tigf) 
— sinh(rig') cosh(Tig') 

cosh(rig) 



exp(TiHi) ^ 



—q^ sinh(r3g^)^ 



cosh(T3g^) J 








exp(r2H2; 



exp(r3F3/) 





V 

^ C0Sh(T2Q'^) 



—q~^ sinh(T2g^) 

V 

^ cosh(T3g^) 



\— sinh(r3g^) 








—q sinh(rig) 
q~^ sinh(rig) cosh(rig) / 



-q^ sinh(T2g^) 

cosh(T2Q'^) — sinh(T2g^) 

cosh(T2g^) 

sinh(T2g^) cosh(T2g^) J 

-g^sinh(T3g^)^ 

cosh(T3g^) — gsinh(T3g^) 

-q~^ sinh(T3g^) cosh(T3g^) 

cosh(r3g3) J 
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And the third group of finite metamorphic operations is 

/e^o \ 
e"^" 



exp(roPo) 



exp(r3P3) 



exp(r3P3 



e- 







\ e^7 

/ cos(r3g^) 



\q~^ sin(r3g^) 



/ cos(r3g'^) 



\q~^ sin(r3g^) 



cos(r3g^j 
sin(r3g^) 









—q sin(r3g^ 
cos(r3g3) 






cos(r3g'^j gsin(r3g'^ 
—q^^ sin(r3g^) cos(r3g'^) 




—q^ sin(r3g'^)^ 



cos(r3g3) / 

—q^ sin(r3g'^)^ 



cos(r3g^) / 



(25) 



(26) 



(27) 



D. Morphological shifting and Jeffrey's third-rank tensor 



Based on the mathematical structure of Ref. [15|, in Ref. |16| a "shifting operation" was 
investigated that changes the radius of a sphere by a given amount R. These operations build 
a one-dimensional Abelian group. The generator of the group, Ti (referred to as G in 



generates the kernel Kr of Ref. 



16|), 



15| upon exponentiation, Kr = exp{RG) = exp(i?Ti). The 



significance of the matrix lies i) in the algebraic structure: K^j ■ K^/ = K^j/ ■ = Kr+^j/, 
and ii) in the fact that it contains the expressions for the four Kierlik-Rosinberg weight 
functions explicitly. 

Jeffrey's third-rank tensor as a central object of Ref. [l6^ is given by the following set of 
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four matrices: 



/l o\ 
10 
10 
\0 1/ 

/o 

-gV(47r) 

1 
\0 1 



\1 

Ti/(87r), T3 = 



Ti 



/o 

1 
Stt 
\0 

-g7(647r2) 





-gV(647r- 







gV(47r) 

1 

\ 

-gV(647r2 





-?V(8vr) 
/ 









(28) 



(29) 



(30) 



where T2 = Ti ■ Ti/(87r), T3 = — g^(Ti)~^/(87r), where (Ti)~^ is the inverse of Ti. 
All commute with each other, [T^, Tj,] = 0. Finite transformations are obtained 
via exponentiation as exp{xu'^u), where Xi^ is the transformation parameter. The fi- 
nite transforms commute with each other and they obey exp(^j^ Xi^Tj,) exp(^j^ xl'^u) — 
exp(Ez. x'l^T^) exp(E,.Xi^T^) = exp(X)^(Xi. + xt)T,.). 



[X,Y]/2 


To 


Ti 


T2 


Ta 


To 














Ti 














T2 














T3 















X- Y 


To 


Ti 


T2 


T3 


To 


To 


Ti 


T2 


T3 


Ti 


Ti 


87rT2 


-|^Ti+T3 


-I^To 


T2 


T2 


-£Ti + T3 


64n'^ 1 I 2 


647r^ ' 1 


T3 


T3 


-I^To 


T. 

647r^ ' 1 


327r2 '0 Stt ' 2 



TABLE IV: Tables for commutator relationships [X, Y]/2 and for products X • Y for the generators 



of transformations. 
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Explicitly, the matrices for finite transformations are given by 

/e'^o \ 
e>^" 



exp(xoTo 



e^» 
\ e^n ) 

I c + qsxi/'i {cq^xi- (ls)/2 



(31) 



exp(xiTi 



/^sxi/(167r) (cg^'xi — 3sg^)/(167r)\ 



exp(x2T2) 



{s + cqxi)/{2q) c-{qsxi)/2 -{3sq + cq\i) / {IGn) 
Ansxi/q Atc{s + cqxi)/q c-{qsxi)/2 
\Att{s - cqxi)/q^ 47rsxi/g (s + cgxi)/(2g) 

[9 + (^?g'X2)/(8vr) -gq\2/{Q^7r^) 

g-gq\2/i8n) 

9X2 g- {gq\2)/{8TT) 



-q^sxi/{16 



TT 



(cg^Xi - sq)/2 
c+ (gsxi)/2 

(32) 

\ 

gq\2/{M7r^) 




V 







9X2 







9 + {9q\2)/{Sn)J 



(33) 



exp(x3T3 



/ C- {q'Sx3)/{lQn) -{STrqS + Cq\3)/{167T) 
S/{2q) - {Cq\,) C + (g35x3)/(167r) 

-qSx,/2 A7iS/q~{Cq\,)/2 

\ 47r^/g3 + CX3/2 -qSx^l2 

g5^X3/(1287r2) -(247rg35 + C q\^) / {12%ti^)\ 

{-2A'KqS + Cg4x3)/(1287r2) g55x3/(1287r2) 

C + (g35x3)/(167r) -(87rg5 + Cq\^)/{lQ'n) 

S/{2q) - Cq^./ilQ-n) C - (g3^X3)/(167r) / 



(34) 



where we have used the short-hand notation s = sin(gxi),c = cos(gxi)) 9 = 
exp(-g2;^2/(8vr)), C = cos(g^X3/(87r)), S' = sin(g3;^3/(87r)). Eq. ([32]) describes Kr when 
setting Xi = R- 

It is an interesting application of the theory outlined in the previous section to try an 
express the Tj, as linear combinations of the metamorphic generators. This can indeed be 
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done with a little algebra, yielding the result: 



To 
Ti 

T2 
T, 



1, 

Fi-Hi 

2 

g^(Po-l) 

Stt 

IGtt 



27r(P3 - + F3 - F3O 3P3 - Ps' - F3 + 3F3, 



2 

P3 + P3' 



-2 + H2 

1287r2 ' 



Fy 3P3 + P3' + F3 + 3F3 

2567r2 



(35) 
(36) 

(37) 
(38) 



III. CONCLUSIONS 



In conclusions we have presented a framework for manipulating four-dimensional vector 
fields u that are defined on an underlying three-dimensional Euclidian space. In real space, 
the relevant operations are application of the Laplace operator and building convolutions. 
These operations turn to multiplication by — and the product operation in Fourier space. 
We have analysed the symmetries that leave the metric ([1]) for the four- vectors invariant. 
This leads to operations that either leave the metric invariant (isometric transforms) or that 
change the metric and hence the morphology that the four-vectors describe (metamorphic 
operations). We have kept the nature of the four- vectors general, i.e. these can taken on 
aribtrary real values. This includes specific geometries (such as spheres considered in Ref. 



161]), but is more general. Whether the transformations presented here help to contruct novel 



DFT approximations is an interesting question for future work. It would also be interesting 
to explore possible connections to the integral geometric framework by Hansen-Goos and 
Mecke Q. 
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Appendix A: Mixed commutator relations of isometric and metamorphic generators 

Here we give further details about the algebra of commutator relations. Tab. |V] gives a 
multiplication table between metamorphic generators as well as anti-commutator relation- 
ships. Tab. I VI I gives products and anti-commutators between mixed pairs of an isometric 
and a morphometric generator. 



X- Y 


Fi F2 F3 


Hi H2 F3' 


Po P3 P3' 


Fi 
F2 
F3 


-qH -F3 q^f2 
-Fa -qH q^fi 
q^\=2 g^Fi qH 


q'^Bo' -P3' -q^B2 
-P3 q^Bo -q^Bi 
-q^D2 -q^Di q^Po 


-Bi -q^D2 q^H2 
-B2 g^Hi -q^Di 
F3, -gSBg, _g6BQ 


Hi 
H2 
Fs' 


-q^Bo> -P3 q^D2 
-P3' -q^Bo q^D^ 
q^B2 q^Bi q^Po 


qH -fy -q^H2 
-F3' qH -g^Hi 
-q^H2 -g^Hi qH 


-Di -q^f2 q^B2 
-D2 q^Bi -q^Fi 
F3 -g^Bo -q'^Bo' 


Po 
P3 
P3' 


Bi B2 F3/ 
q^D2 q^Hi q^Bo, 
q^H2 q^Di q'^Bo 


Di D2 F3 
-q^f2 -g'Bi q^Bo 
-q^B2 -q^fi <7'Bo' 


1 P3' P3 

P3, _g6p^ 

P3 -g'Po -g'l 



{X,Y}/2 


Fi 


F2 


F3 


Hi 


H2 


F3' 


Po 


P3 


P3' 


Fi 




-F3 


q'f2 





-P3' 











q'^2 


F2 


-F3 


-q'l 


q'fl 


-P3 











g^Hi 





F3 


<Z^F2 


q'fl 


q^l 








g'^Po 


F3' 








Hi 





-P3 





q'l 


-F3' 


-q'^2 





-q'f2 





H2 


-P3, 








-F3' 


q'l 


-g'Hi 








-q'fi 


F3' 








q^Po 


-q'y^2 




q^l 


F3 








Po 








F3' 








F3 


1 


P3' 


P3 


P3 





g^Hi 





—q r2 








P3' 


-q^l 


-q'Po 


P3' 


q'^2 











-'Z^Fi 





P3 


-q^Po 


-q^l 



TABLE V: Top: Multiplication table X-Y for products of the generators of metamorphic operations. 
Bottom: Anti-commutator relations {X, Y}/2 = (X • Y + Y • X)/2 for generators of metamorphic 
transformations. X denotes a matrix of the leftmost column, Y one of the top row. 
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X • Y 


Fi 


F2 F3 


Hi 


H2 


F3 


Po 


P3 


P3 


Bo 


Di 


H2 P3 


Bi 


F2 


P3 


B'o 


F^ 


F3 


B2 


-f:3 


-P09' Big4 


-P'a 




-Fig4 


-F2 


Dig^ 


-Hig^ 


D2 


-P3 


-%q^ Hig4 


-F3 






-H2 


Fig^ 


-Bi<?^ 


B^, 


Hi 


D2 P3 


Fi 


B2 


P'3 


Bo 


F3 


F^ 


Bi 




-F'3 B2g2 


Bog^ 


-P3 


-F2g2 


-Fi 


-H2g2 


□29^ 


Di 


-Bo?^ 


-P3 H2g2 


Pog' 


-F3 


-D2g^ 


-Hi 


-B2g^ 


F29^ 



YX 


Fi 


F2 


F3 


Hi 


H2 


F'3 


Po 


P3 


P' 

^3 


Bo 


Di 


-H2 


— P' 


Bi 


-F2 


-P3 


B^, 


-F^ 


-F3 


B2 


F's 


Pog^ 


Big4 


P3 


%q' 


Fig^ 


F2 


Dig^ 


Hig^ 


D2 


P3 


-Kq' 


-Hig4 


F3 


-Pog' 




H2 


-Fig^ 


-Big^ 


B'o 


-Hi 


D2 


-P3 


-Fi 


B2 


-P^ 


Bo 


-F3 


-F^ 


Bi 


P09' 


F^ 


B2g^ 


Bog^ 


P3 


F2g' 


Fi 


H2g^ 




Di 


-Bog^ 


P^ 


-H2g2 


-Pog' 


F3 


-D2g^ 


Hi 


-B2g^ 


-F2g^ 



{X.Y}/2 


Fi 


Fo 


F:^ 


Hi 


H2 


F'5 


Pn 


p.i 




Bo 


Di 








Bi 








B^ 








B2 








Big^ 





Kq' 








Dig^ 





D2 





-B'o?^ 











-Dig^ 








-Big4 







D2 








B2 





Bo 








Bi 








B2g^ 


Bo?^ 














□29^ 


Di 


-Bog^ 














-D2<z' 





-B2g^ 






TABLE VI: Top: Multiplication table between generators of isomorphisms, X, and generators of 
metamorphisms, Y. Middle: Reverse product order, Y-X. Bottom: Anti-commutator relationships 
for the same pairs. 
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